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Abstract 

This article is meant as a mathematical appendix or comment on pp. We first con- 
sider the notion of transcritical bifurcations of fixed points of general area-preserving 
maps, and then adress some questions related to [1] on bifurcation in Poincare maps 
of 2-dimensional Hamiltonian systems. 



1 Rank-l-bifurcations 

Let different iable (C°°) functions Q = Q(q,p,e) and P = P(q,p,s) be defined on 
an open neighborhood of the origin in (q,p, e)-space M 3 . We assume (P,Q) to be 
symplectic in the (g,p)-coordinates, which means that 

det f Q q (q,p,e) Q P (q,p,e) 

where we have written the partial derivatives as |p =: Q q etc. We then speak of 
(Q, P) as describing a symplectic family. Our bifurcating fixed point shall be the 
origin of the (g,p)-plane at the value e = of the parameter, so Q(0, 0, 0) = and 
P(0,0,0) = 0. At the bifurcation point (0,0,0), we assume the Jacobian matrix 
with respect to the variables (q,p) to have the eigenvalue +1. Otherwise, by the 
implicit function theorem, the fixed point set could be parametrized locally by e and 
hence would not bifurcate in the way we wish to study. We will exclude, however, 
the exceptional case that the Jacobian matrix equals the identity matrix. So by 
assumption, the eigenspace is a 1-dimensional subspace of the (g,p)-plane. We can 
always adjust the canonical coordinates, for example by a simple rotation, to have 
the eigenspace as the g-axis. Then at (0, 0, 0) we have 



Qq Qp \ f 1 Qp 



Pq Pp I V01 



and Q p ^ (2) 



in these coordinates. The total fixed point set 

F ■= {(q,p,e) | Q(q,p,e) = q\ P(q,P,e) = p] (3) 

is the inverse image of the origin (0, 0) in M. 2 under the map (Q — q,P — p), and the 
Jacobian matrix of this map at (0, 0, 0) is 



Qq-l Qp Qe\ _ ( Q p Q 



P q P p -1 P e \ P e 



(4) 



If we were to look for generic bifurcations, we could now impose the condition that 
this matrix be of rank 2, or P £ ^ at (0, 0, 0). Then F would be locally a smooth 
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1- dimensional submanifold, tangent to the g-axis at the bifurcation point, and if we 
require, as a further generic condition, e\F to have a nondegenerate extremum at 
this point, we would arrive at the notion of an extremal fixed point or saddle-node 
bifurcation in the sense of K. R. Meyer [2]. 

Instead, in this note we will be interested in bifurcations with matrix (4) being 
of rank 1. Among these rank- 1 -bifurcations we will distinguish various types 
by further conditions. For a convenient coordinate description, let us adjust the 
canonical (q, p)-coordinates one step further by an e-dependent translation 

1= q (5) 
p = p — ce 

where the constant c is defined by c := —Q e /Q P at (0, 0, 0). Then the g-axis still is 
the eigenspace at the bifurcation point, and the Jacobian matrix of (Q — q, P — p) 
at (0, 0, 0) is simplified to 



Qq 1 Qp Qe 

p~ P 1 p 

1 q 1 p J- 1 e 



)- : 


Qp 










o / 



(6) 



So let (q,p) be such coordinates to begin with. To have a name for them, we 
introduce the following terminology. 

Definition 1: Symplectic coordinates (q,p) shall be called adapted coordinates 
for a rank-l-bifurcation at (0, 0, 0), if 



Qq-i Qp Qe\ = ( Qp 

P q Pp-l Pe I 



(7) 



at the bifurcation point, and Q p ^ there. □ 

Then the inverse image of under Q — q, let's denote it by 

X := {(q,p,e) | Q(q,p,e) = q}, (8) 

will locally at = (0,0,0) be a smooth surface, with the (q, e)-plane as tangent 
plane T X at this point. The fixed point set F is contained in X, it is the inverse 
image 

F = {(q,p,e)eX\P(q,p,e)=p}, (9) 

of under the restriction (P — p)\X of P — p to X. As (0, 0, 0) is a critical point of 
this restriction, we naturally turn to its Hessian quadratic form on the (q, e)-plane 
for information about the local behavior of P — p on X. And since (0, 0, 0) is critical 
not only for the restriction, but also for the function P — p itself, we know that this 
quadratic form is simply described by the Hessian matrix of P — p, hence of P, at 
(0, 0, 0) with respect to the (q, e)-coordinates. This is how this matrix enters the 
following definition. 

Definition 2: A rank-l-bifurcation shall be called regular, if in adapted coordi- 
nates the Hessian 

p qq p Q£ ) at (0,0,0) (10) 
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of P with respect to q and e is nondegenerate. Depending on whether it is definite 
or indefinite, we speak of the rank-l-bifurcation as being definite or indefinite. 
An indefinite rank-l-bifurcation will be called a cross-bifurcation, and a cross- 
bifurcation is called trans critical, if the g-axis is not contained in the zero set of 
the quadratic form, that is if P qq ^ at (0, 0, 0), in adapted coordinates. □ 

The definition does not depend on the choice of the adapted coordinates. As can 
be shown, the Hessians from different choices of adapted coordinates are equivalent 
(up to sign) as quadratic forms. 

2 The fixed point set: branches and traces 

Let a regular rank-l-bifurcation in adapted coordinates be given. We now denote 
the restriction (P — p)\X by ip : X — > R. So we are interested in ^ _1 (0) C X, 
because this is the fixed point set F of the bifurcation. The Hessian quadratic form 

Hess^, : T X -> R (11) 

of ip at the bifurcation point is given by the matrix (10). In coordinates (q,e) on 
X, Taylor expansion to second order of ip at this point gives ip ~ |Hess^ up to 
higher order terms. But since the Hessian is assumed to be nondegenerate, we can 
do better than that. By the Morse Lemma, see for instance [3j, we can find local 
coordinates on X, approximating (q, e) in first order, in which ip actually coincides 
with |Hess^. More precisely, there is a diffeomorphism / : Q — > Q' from an open 
neighborhood Q of the origin in the tangent plane of X to an open neighborhood Q' 
of the bifurcation point in X itself, such that /(0) = and 

iP(f(v)) = §Hess^) (12) 

for all v G O, and the differential dfo '■ TqX — > TqX is the identity. Thus up to this 
diffeomorphism, the fixed point set of the bifurcation looks locally the same as the 
zero set Hess^ 1 (0) of the Hessian, which is a single point in the definite and a pair 
of straight lines in the indefinite case. Summing up: 

Proposition 1: In a sufficiently small neighborhood of a regular rank-l-bifurcation 
point (0,0,0) in the (q,p,e) -space, the fixed point set consists of a single point if 
the bifurcation is dehnite and it is the union AU B of two smooth 1-dimensional 
submanifolds, intersecting at the bifurcation point, if the bifurcation is indefinite 
that is in the case of a cross-bifurcation. Moreover, if a and b denote the tangents 
to A and B at the bifurcation point, then a ^ b, and the intersection of the plane 
spanned by a and b in the (q,p,e)-space with the plane defined by fixing the bi- 
furcation parameter value that is with the (q,p) -plane e = is the 1-dimensional 
eigenspace at the bifurcating fixed point. The bifurcation is transcritical, if and only 
if this eigenspace is different from a and from b. □ 

Recall that in adapted coordinates the eigenspace is the g-axis and the plane spanned 
by a and b is T X, the (q, e)-plane. 

Now let for a moment (Q(p,q,e), P(p,q,e)) denote any local symplectic family 
with a fixed point at (0,0,0), without assuming it to be rank-l-bifurcating or bi- 
furcating at all. If there is a different iable map a : / — > R 2 , defined on an open 
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interval i" C R around 0, written in coordinates as a(e) = (q(e),p(e)), such that 
a(0) = (0,0) and 

= 9(e) and P(g(e),p(e), e) = p(e) (13) 

for all £ £ I, then we call its graph 

A:={(q(e),p(e),e)\eeI}cR 3 (14) 

a fixed point branch of (0, 0, 0), regardless of what other fixed points of the sym- 
plectic family might exist. For instance, if a transcritical bifurcation is considered in 
a suitable small neighborhood, then its fixed point set will consist of two branches A 
and B as in Proposition 1, but if an indefinite rank-l-bifurcation is not transcritical, 
then only one of its fixed point lines, say A, will be a branch, while B, being tangent 
to e = 0, cannot be parametrized differentiably by e. 

If again A is a fixed point branch of (0, 0, 0) in some symplectic family, then the 
trace of the Jacobian matrix, taken at each point of the branch, defines a differen- 
tiable real valued function Tr^ : I — > R, so 

Tr A (e) :=Q,(g(e),p(e),e) + P p (g( e ),p( e ),e). (15) 

Then the eigenvalue at (0,0,0) is 1 if and only if Tr^(0) = 2. Under what addi- 
tional condition will (0, 0, 0) be an indefinite rank-l-bifurcation point, i.e. a cross- 
bifurcation? 

Proposition 2: If a fixed point branch A of (0, 0, 0) in a symplectic family satisfies 
Tr^(0) = 2, then (0, 0, 0) will he a cross-bifurcation point if and only if Tr^(0) ^ 0. 

Proof: So let a fixed point branch A of (0, 0, 0) in a symplectic family (Q, P) be 
given, and Tr A (0) = 2 be assumed. We have to show (a): If Tr^(0) ^ 0, then (0, 0, 0) 
is a cross-bifurcation point, and conversely (b): If (0, 0, 0) is a cross-bifurcation point, 
then Tr^(0) ^ 0. 

Proof of (a): We first choose the symplectic (q, p)-coordinates in such a way, that 
the g-axis is contained in the eigenspace, and the e-axis q = p = is tangent to the 
branch at (0,0,0). Then we already have 

( Q q -i Q P Qe \ _ ( Q p \ 
V P q P p -i P £ ) V J ' 

at (0, 0, 0), but we do not yet know if Q p ^ 0, and we know nothing about the Hessian 
of P, except P ee (0, 0, 0) = 0, which follows from differentiating P(q(e),p(e),e) = p(e) 
twice at e — 0, note g'(0) = p'(0) = by our choice of coordinates. So it remains 
to show that Q p ^ and P q£ ^ at the bifurcation point. Now since (Q, P) is a 
symplectic family, we have 

det ( «■ «;)=!, (17) 

everywhere, not only at fixed points. So if u denotes any of our three variables q, p 
or e, we obtain 

det ( Qqu Q p p \+ det ( Qq Qpu \=0 (18) 
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everywhere, and thus at the special point (0, 0, 0) we get 

det(« : ?') + d*(J «-) 

For later reference, let us call this the 

Determinant derivative formula: If Q q = P p = 1 and P g = at some point in 
a symplectic family , then if u is any of the three variables q, p, or e we have 

qu Ppu Qp-Pqu 

at this particular point. 

Presently we apply the formula for u := e. Since Tr^(0) = Q q£ (0, 0, 0) + P pe (0, 0, 0) 
because of q'(0) = p'(0) = 0, we have Tr^(0) = Q p P q£ at (0,0,0), so Q p P q£ ^ by 
our assumption Tr^(0) ^ 0, and (a) is proved. 

Proof of (b): Now we assume P qq P ££ — P q£ < at (0,0,0) in adapted coordinates, 
and we have to show that this implies Tr A (0) ^ 0. The branch A must be tangent 
to the (g,e)-plane at the bifurcation point, which means p'(0) = 0, and a change to 
adapted coordinates of the type 

q = q-c x e 
p = p 

with a suitable constant c\ will even make the e-axis tangent to A, with the Hessian 
still being indefinite in the new coordinates. The trace function Tr^(e) remains 
unchanged anyway, as it is independent of the coordinate choice. So we may assume 
g'(0) = p'(0) = for the branch from the start. As before, this implies Tr^(0) = 
Q q£ + P p£ and P ££ = at (0,0,0). Since we have a rank-l-bifurcation, we know 
Q p 7^ 0, and since it is regular, P q£ ^ follows from P ££ = at (0,0,0). The 
determinant derivative formula is applicable and shows Tr^(0) = Q p P q£ and hence 
Tr^(0) 7^ 0, which completes the proof of (b) and of Proposition 2. □ 

As we have seen in Proposition 1, locally at a transcritical bifurcation point the 
fixed point set consists of two branches A and B. Since the functions Tr^e) —2 and 
Trg(^) — 2 both change sign at e — 0, both branches change their stability properties 
from elliptic to hyperbolic or vice versa, and in fact in opposite directions. More 
precisely 

Proposition 3: If A and B are the branches of a transcritical bifurcation point, 
then Tr' A (0) + Tr' B (0) = 0. 

PROOF: In adapted coordinates we have 7/4(0) = p#(0) = anyway, since locally 
the fixed point set is contained in the surface X tangent to p — 0. Changing the 
adapted coordinates by a suitable transformation of the type (21), we can also obtain 
<?a(0) +<?b(0) = 0- I n these coordinates, we have P q£ (0, 0, 0) = 0. On the other hand, 
Tr / A (0)+Tr' B (0) = g w ^(0)+Q, e +P M ^(0)+P pe +g M g^(0)+Q, e +P M ^ B (0)+P pe = 
2(Q (?e + Ppe), but this is 2Q p P q£ by the determinant derivative formula, and the 
proposition follows. □ 



= 0. 



(19) 



(20) 



(21) 
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3 Fork-like bifurcations 

Let us now have a look at those cross-bifurcations, which are not transcritical. Then 
only one of the two fixed point lines A and B of Proposition 1 will be a branch in 
the technical sense (14), while the other, say B, is tangent to the g-axis. Thus the 
function B — > R given by the restriction e\B of the family parameter to B is critical 
at the bifurcation point. It may in fact happen that e\B = 0, the most degenerate 
possibility. But here we introduce a terminology for the least degenerate case. 

Definition 3: Let a cross-bifurcation which is not transcritical be given. Then if 
the restriction e\B of the bifurcation parameter to one of the local fixed point lines 
has a nondegenerate extremum at the bifurcation point, we speak of a fork-like 
bifurcation. □ 

It is the same condition on B as is required of the fixed point set in the definition 
of an ordinary saddle-node rank-2-bifurcation. But in our case, the additional fixed 
point branch A is present, while in a saddle-node bifurcation the fixed point set is 
locally just a single line. 

Since B in the non-transcritical case is tangent to the eigenspace, it can be para- 
metrized as (<?, P_b(<?), £b(<?)) by q in adapted coordinates, with p' B (0) = £ B (0) = 0. 
The fork-like condition then just says £ B (0) ^ 0. In the surface X of (8), the line B 
locally looks like a parabola up to higher order terms, tangent to e — 0. 

The name is meant to refer to the well-known pitchfork bifurcations in symplectic 
families. These are bifurcations of fixed points that have eigenvalue —1 and as 
such do not fall under our heading. But if we consider the iterated family instead, 
(Q(Q(q,p, s), P(q,p, £),£), P(Q(q,P, s), P(q,P, £),£), z) in coordinates, then the fixed 
point gets the eigenvalue +1 and the iterated family bifurcates it fork-like. But being 
iterations, these are very special fork-like bifurcations. For instance does the original 
family provide a sort of symmetry, which a fork-like bifurcation in general will not 
have. 

The trace along the fixed point line B of a fork-like bifurcation is of course well- 
defined as a function Tr B on B itself, but not like Tr A as a function of e, since for 
small e^Owe find either none or two points on B, and these two may in fact have 
different traces. But in an natural way, Tr' B (0) is still defined and in a fixed relation 
to Tr^(0), as the next proposition shows. 

Proposition 4: Let A be the fixed point branch of a fork-like bifurcation and B 
the fixed point line tangent to e = 0. Then there exists the limit 



lim 

q^0 



Tr B (q,p B (q),e B (q))-2 
e B {q) 



: Tr' B (0) 



(22) 



and it satishes 



Tr' A (0) + 



iTr' B (0) = 0. 



(23) 



PROOF: To determine the limit, we will use second order Taylor expansion of 
Q q (q,p,e) + P p (q,p,e) in adapted coordinates (7). All partial derivatives now taken 
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at (0, 0, 0). By the determinant derivative formula we know 

{Qq "I" Pp)q = QpPqq 

{Qq ~\~ Pp)p = QpPqp (24) 
{Qq Pp/E QpPqe j 

and P qq = since the bifurcation is not transcritical. Moreover, 

Tr^(0) = (Q q + P P ) q q' A {0) + {Qq + P p )pP'a{0) + {Qq + P P )e (25) 
and p' A {0) = 0, so from (24) and P qq = we also get 

Tr^(0) =Q p P q£ . (26) 
Thus we know the linear terms, and the quotient becomes 



(27) 



Tr B (q,p B (q),e B (q)) -2 p B (q) 

H — • (higher order terms) 

But not many of the higher order terms will contribute to the limit, since 

i im ^44=4^ and ^^ = 2, (28) 

i^e B {q) 4(0) e B (q) 4(0)' 1 ; 

and as an intermediate result we obtain that the limit Tr^(0) exists and that 

Tr' B (0) = Tr' A (0) + 3 P ^Pb(°) + Qm+Jjw (2Q) 

£ b(0) 

To prove (23), we have to prove that the quotient on the right hand side of this last 
equation (29) equals — 3Tr^(0), and thus by (26) what remains to be shown is 

3Q p P q£ e B {0) + Q p P qp p" B {^) + Q qqq + P pqq = . (30) 

Now we use the fixed point property 

Q{q,PB{q),e B (q)) = q 
P{q,PB{q),e B (q)) =p B {q) 

of B. Differentiating the first equation twice and the second three times at q — 0, 
we get 

Qqq + QpP' B {0)=0 (32) 

and 

Pqqq + 3P q£ 4(0) + 3P qpP ' B (0) = 0. (33) 

Inserting this into equation (30), which we have to prove, we find that (30) is 
equivalent to 

Qqqq ~t~ Ppqq QpPqqq 2Q qq P pq (34) 

at (0,0,0). But this turns out to be a consequence of a second order determinant 
derivative formula: we start from the symplectic property (17), apply d 2 /dq 2 , put 
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in what we know about the first and second order partial derivatives of Q and P at 
(0, 0, 0), and out comes (34). Thus Proposition 4 is proved. □ 

From the calculations of the proof, let us preserve the formula 

4 (0) = gggjkzgejk ? (35) 

which is a consequence of (32) and (33). Since in its derivation the non- vanishing 
of e' B (0) was not used, we have as a 

Corollary: A non-transcritical cross-bifurcation is fork-like if and only if 

3QqqPqp 7^ QpPqqq (36) 

at (0, 0, 0) in adapted coordinates. □ 
4 Fixed point branches given by librating orbits 

Now let H : M — > K be an autonomous Hamiltonian on a 4-dimensional symplectic 
manifold and 7 : R — > M a periodic orbit with period T > of the Hamiltonian 
flow. Choose any 3-dimensional submanifold E C M which is being intersected 
transversally by 7 at time i = 0. Then on a sufficiently small neighborhood U C E 
of 7(0) in E, the Poincare map Poinc : £7 — > E is well-defined by following the orbits 
until they hit E again after travelling approximately the period time T of 7. Let 
£"o := if (7(0)) be the energy of the fixed point 7(0) and write Ue '■= UPiH~ l (E) and 
Eg := En H~ l (E). After making E smaller if necessary, the individual Ue and Eg 
will be nondegenerate subsurfaces of M, with the Poincare map defining a symplectic 
map Ue — > E# for each i?. Now we introduce a coordinate £ for the energy, like 
e = E — Eq, and extend it to local coordinates (q,p,e) for E such that (q,p) are 
symplectic coordinates on each E#. Then the Poincare map will be a symplectic 
family, described by two functions Q(q,p,e) and P(q,p,e) in these coordinates. All 
this is of course well-known, recalled here only to introduce notation. 

Often used is the following sort of 'automatic' choice of E and the coordinates q 
and p. Let H = H(x, y,p x ,P y ) be given on M = IR 4 . If the periodic orbit 7 satisfies 
|)(0) 7^ 0, then the 3-dimensional subspace y = yo := y(0) in M 4 can be taken for a 
start to find E. Of course, far away from the point 7(0) this space may have bad 
properties with respect to the Hamiltonian flow. But if we choose E as a sufficiently 
small open neighborhood of 7(0) in this 3-space, then not only the Poincare map 
and the energy surfaces E^ will be defined as described, but also the projection of 
each Tje to the (x,pa;)-plane will be symplectic and injective, which means that x 
and p x can be used as the symplectic coordinates q and p on E E . 

The transcritical bifurcations studied in [T] are related to straight-line librating 
orbits of Hamiltonians of the form 

H{x,y,p x , Py ) = \pl + \p 2 y + V{x,y). (37) 

Let 7 be such a straight-line periodic orbit, without loss of generality projecting 
to the y-axis, oscillating there between two values y± < y%, with x(0) = and 
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yi < y(0) < %)2 and an energy E . As just recalled, we use 



q := x 



Px (3£ 
E — Eq 



as coordinates to describe the 'automatic' Poincare map by a symplectic family 
(Q(q,p,e), P(q,p,s)), with the fixed point (0,0,0) corresponding to the orbit 7. 
What do we know about this symplectic family? 

Straight-line librating orbits in Hamiltonian systems of type (37) always come 
in one-parameter families. The given orbit 7 satisfies 

dV 

p x (t) = - — (0,y(t)) = 0, (39) 

hence at least we know ^r(0, y) = for all y with y%<y < yi- If V is polynomial 
or real-analytic, this implies 

dV 

— {0,y) = for all yeR. (40) 
ox 

In this case the whole (y,p y ) -plane x = p x = in M 4 is invariant under the Hamilto- 
nian flow, and the flow lines are just those of the 1-dimensional Hamiltonian system 
on the (y, p 2/ )-plane given by 

h(y, Pv ) = ffi + V{0,y). (41) 

By standard regularity arguments, the orbit 7 must be embedded in a family of 
neighboring closed orbits of this 1-dimensional system, one for each energy in an 
interval (Ei,Ez) with E\ < Eq < E 2 . Thus we have a family of librating orbits 
of the original 2-dimensional system on the y-axis ovei0 (E 1 ,E 2 ). In our Poincare 
symplectic family it corresponds to a fixed point branch 

A := {(0,0, e) | e 1 < e < e 2 }. (42) 

For a branch A, the trace function Tr^ : (£1,62) ^ K is defined, and if it happens 
to be that Tr^(0) = 2 and Tr^(0) 7^ 0, then by Proposition 2 we have a cross- 
bifurcation point, which may or may not be transcritical. Examples of both cases 
have been presented numerically in pQ. 

Hamiltonians of type (37) have time-reversal symmetry, and in jl] the destruction 
of a transcritical bifurcation by a symmetry breaking perturbation is described. 
What role does time-reversal symmetry play in transcritical bifurcations? Could we 
preserve a transcritical bifurcation under a symmetry breaking perturbation, and 



Here we assumed (40). What if §£(0,y) may be non-zero for y > yi or y < yi? Then the 
neighboring closed orbits on the inside of 7, those with energies E\ < E < E$, still are librating 
orbits of the original system as well, and define a fixed point branch in the Poincare family with 
q = p = and £1 < e < 0. The orbits of the 1-dimensional system which are on the outside of 7, 
with energies Eq < E < E%, need not be orbits of the 2-dimensional system. But then we simply 
shift attention from Eq to an energy Eq between E\ and Eq. So in any case, we may assume a 
fixed point branch A as in (42) in the Poincare symplectic family be given, describing a family of 
neighboring straight-line librating orbits on the y-axis. 
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conversely, can a transcritical bifurcation be destroyed by a symmetric perturbation? 
The answer to the first question is yes, as we are going to show now. 



5 Transcritical bifurcation in unsymmetric Hamiltonian systems 

Again we start from a Hamiltonian H(x,y,p x ,p y ) = \p\ + \p^ + V(x,y) with a 
straight-line librating orbit 7, embedded in a family of such orbits, all projecting 
to the |/-axis and thus constituting a fixed point branch A = x x (£1,62) in the 
Poincare symplectic family as described above, with (0, 0, 0) representing 7. For 
convenience and with little loss of generality we assume (40), that was ^(0, y) = 
for ally G R. But our main assumption now will be that (0, 0, 0) is a cross-bifurcation 
point in the Poincare family. We then speak of 7 as of a cross-libration or of a 
transcritical libration, if this cross-bifurcation happens to be transcritical. 

Now let the Hamiltonian depend on an additional small parameter S. Such a 
function H(x,y,p x ,p y ,S) will be called a perturbation of H or of the transcriti- 
cal libration, if H(x,y,p x ,p y ,0) = H(x,y,p x ,p y ). Then we have a two-parameter 
Poincare symplectic family 

(Q,P) = (Q(q,p,e,S),P(q,p,e,S)), (43) 

defined on an open neighborhood of (0, 0, 0, 0) in the (q, p, e, 5)-space R 4 in the usual 
way, with q and p from x and p x . 

Definition 4: We say that a perturbation if of a cross-libration ist cross-preser- 
ving, if there are differentiable functions q(S), p(5) and s(S), defined on an interval 
(— So, S ), with q(0) = p(0) = e(0) = 0, such that for any fixed S G (—So, So) the 
point (q(S),p(S),e(S)) is a cross-bifurcation point in the symplectic e-family given 
by (Q, P) at the fixed S. □ 

Note that a transcritical libration will stay transcritical under a cross-preserving 
deformation, for small enough S - 

There is a simple strategy to find cross-preserving perturbations. All we have to 
do is to make sure that for all sufficiently small S the (y , p y )-pl&ne OxRxOxR 
is invariant under the Hamiltonian flow on R 4 that is defined by H for fixed S. In 
other words, if 

OH OH 

— (0,y,0,p y ,S) = and — (0, y, 0,p y , S) = (44) 

for sufficiently small S and all (y,p y ), then H will preserve the transcriticality. Why? 
Because then the closed orbits neighboring 7 will define a 2-parameter fixed point 
branch A of (Q, P) given by functions q(e, S) and p(e, S), defined on a neighborhood 
of (0, 0) in the (e, <5)-plane. We know 

OTr— 

Tr^(0,0) = 2 and -^(0,0)^0 (45) 

by assumption, thus by the implicit function theorem, we get a differentiable function 
e(S) with e (0) = and 

OTr— 

Tt^(s(S),S) = 2 and -^( £ (S),S)^0 (46) 
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for sufficiently small 5. By Proposition 2, each (q(e(S),S),p(e(S),S),e(S)) is a cross- 
bifurcation point in the symplectic family of the corresponding fixed 5 (and if it is 
transcritical for 5 = 0, it must remain transcritical for sufficiently small 5). So the 
strategy is sound, and we use it to derive the following proposition. 

Proposition 5: Let H(x,y,p x ,p y ) = \p\ + \p 2 y + V(x,y) be a Hamiltonian with a 
cross-libration on the y-axis and dV/dx(0, y) = for all y e R. Let F(x,y,p x ,p y ) 
be an arbitrary differentiable function with 

OF OF 

— (0,y,0,p y ) = and — (0, y, 0,p y ) = (47) 

Then 

H(x, y, p x , p y , 5) := H(x, y, p x , p y ) + 5 F(x, y, p x , p y ) (48) 
is a cross-preserving perturbation, since it obviously satisfies (44). □ 

It is a matter of definition, which Hamiltonians should be counted as 'unsymmetric' 
in this context. Making H unsymmetric by a bump far away from the transcritical 
libration would not illuminate the relation between symmetry and transcriticality. 
But Proposition 5 shows transcritical bifurcation to occur in the Hamiltonian system 
H(x,y,p x ,p y ) + 5 F(x,y,p x ,p y ) on M 4 for fixed sufficiently small 8^0. Since the 
only condition on F(x,y,p x ,p y ) is the vanishing of the partial derivatives by x and 
p x along the (y, p J/ )-plane, we may fairly say yes, transcritical bifurcation can occur 
in unsymmetric autonomous 2-dimensional Hamiltonian systems. 



6 Destruction of cross-bifurcations 

Once more we start with a cross- or a transcritical libration 7 on the y-axis of a 
Hamiltonian system H(x,y,p x ,p y ) = \p\ + \p 2 y + V(x, y) with ^(0,y) = for all 
y G R. This time we ask by which perturbations of H the cross-property can be 
destroyed. By destruction we mean more than just non-preservation: 

Definition 5: We say that a perturbation if of a cross-libration is cross- destroy- 
ing, if there is a neighborhood of (0,0,0,0) in (q, p, e, <5)-space, in which (0,0,0,0) 
is the only cross e-bifurcation point of (Q, P). □ 

We will first prove a general destruction criterion for cross-bifurcation in symplectic 
families. 

Proposition 6: Let (Q(p,q,e), P(p,q,e)) be a symplectic family with a cross- 
bifurcation at (0,0,0). Extend it to a symplectic family (Q(p,q,e,5), P{p,q,e,5)) 
depending on a second parameter S, with Q(p, q, e, 0) = Q(p, q, e) and P(p, q, e, 0) = 
P(p,q,e). Then if the vector (Q s (0, 0, 0, 0), P s (0, 0, 0, 0)) does not belong to the 
eigenspace of (Q, P) at the bifurcation point, there is a neighborhood of (0, 0, 0, 0) 
in the (q,p,e,S)-space, in which it is the only e -cross-bifurcation point. 

Proof: The assumption simply means that at (0, 0, 0, 0) the Jacobian matrix 

/ Qq — 1 Qp Qe Q_6 \ (aq\ 
\ P q P p -1 P £ PS J 1 } 
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of (Q — q, P — p) is of rank 2, as can be seen most easily in adapted coordinates, 
where the condition reduces to Ps ^ 0. So the fixed point set F, locally at (0, 0, 0, 0), 
is a smooth surface, and since it contains the two fixed point lines A and B of the 
given cross-bifurcation, its tangent plane at the point must be spanned by a and b 
und thus be the (q, e)-plane in adapted coordinates. In particular, the restriction 
S\F is singular at (0, 0, 0, 0). 

To prove the proposition, it would be sufficient to show that this is an isolated 
singularity of S\F. Because then at any other fixed point (qo,Po, £o, $o) i n a neigh- 
borhood, the set of fixed points with 8 — So will be locally a smooth 1-dimensional 
submanifold of F and hence the point can't be an e-cross-bifurcation point. 

To find out, we calculate the matrix of the Hessian quadratic form of 8\F at the 
singularity. Straightforward calculation using Lagrangian multipliers gives 

Hess 5F (0,0,0,0) = -=U (50) 

for this matrix, which is now seen to be nondegenerate by assumption, hence the 
singularity is isolated and Proposition 6 is proved. □ 

Let 7 be a cross-libration on the y-axis of a Hamiltonian system 

H(x, y, p x , p y ) = ^p 2 x + ^p 2 y + V (x, y) (51) 

with ^(0, y) = for all y G R. For an attempt to apply Proposition 6 to pertur- 
bations 

H(x, y, p x , p y , 5) := H(x, y, p x , p y ) + 5 F(x, y, p x , p y ) (52) 

with an arbitrary perturbation term F(x,y,p x ,p y ), we will have to understand how 
(Q s (0, 0, 0, 0), P s (0, 0, 0, 0)) does depend on F(x, y,p x ,p y ). 

Proposition 7: Let T > denote the period of the cross-libration 7 and y(t) its 
y-component. Define the function f(y) by 

f(y)--=^(o,y). (53) 

and let gi(y,p y ) and g 2 (y,P y ) be given by 

dF 

9i(y,p y ) ■= T—(0,y,0,p y ), 

opx / 54 x 

dF 

92(y,Py) ■= —^(0,y,0,p y ). 
Let (£(t),rj(t)) be the solution of 

i -V = 9i(y(t),y(t)) ( „\ 

v +f(y(tm = 92(y(t),y(t)) 

with the initial condition £(0) = rj(0) = 0. Then H satisfies the destruction criterion 
of Proposition 6 if and only if (£(T), i](T)) is not contained in the eigenspace of the 
Poincare map of the undisturbed system at the bifurcation point. □ 
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What would it take to apply Proposition 7 numerically? First of all, one has to know 
the component y(t) and its derivative y(t) = p y {t) of the cross-librating orbit 7. 
Secondly we will need the fundamental system ((p(t),i/j(t)) of the homogeneous 
equation 

£+ /(?/(*))£ = (56) 

with initial conditions <p(0) = ^(O) = 1 and 0(0) = ^(0) = 0, and finally we would 
have to calculate the two numbers 



Cl (T) := J(^(t)gMt),y(t))-^(t)g 2 (y(t),y(t))dt 
C2(T) :- 





T 



(57) 



f(tp(t) gi (y(t),y(t)) - <p(t)g 2 (y(t),y(t)) dt. 




This then is all we need, because as can be shown, the Jacobian matrix of the 
Poincare map at the bifurcation point turns out to be 



Qq Qp 



<p(t) vco 

<p(T) ip{T) 

and standard calculation of (£(T),rj(T)) leads to the following 



P P 



(5f 



Corollary: The perturbation term F(x,y,p x ,p y ) satisfies the destruction criterion 
of Proposition 6 if and only if 



<p(T) - 1 i>{T) 
<p(T) 4>(T) - 1 



ci(T) 
c 2 (T) 



(59) 
□ 



Application to F(x,y,p x ,p y ) = xp y — yp x leads to the condition 



ip{T) - 1 v>CO 

<p(T) ip(T) - 1 



V 



jy(t)ip(t)dt ^ 



fy(t)<p(t)dt 




(60) 



For a more systematic approach to the calculation of partial derivatives of the 
Poincare map in connection with straight-line librations, see [3]. 
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